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FOR EWORD 


(U)  This,  the  third  Quarterly  Technical  Report  under  Contract  No.  F 04611- 

67-C-0034,  covers  the  work  performed  from  1  June  through  31  August  1967. 
This  contract  with  the  Huntsville  Division  of  Thiokol  Chemical  Corporation 
was  initiated  under  Air  Force  Rocket  Propulsion  Laboratory,  Research  and 
Technology  Division  Project  Number  3148.  It  is  being  accomplished  under 
the  technical  direction  of  R.  W.  Bargmcver,  1/Lt.  ,  USAF  of  the  Research 
and  Technology  Division,  Air  Force  Systems  Command,  United  States  Air 
Force,  Edwards  Air  Force  Laso,  California  93523. 

(U)  Dr.  David  A.  Flanigan  of  Thiokol's  Research  and  Development  Depart¬ 
ment  is  the  Principal  Investigator  and  Mr.  Carl  J.  Whelchel  of  the  Project 
Management  Directorate  is  Associate  Program  Manager  for  this  program. 

Full  authority  for  the  management  control  of  this  program  is  the  responsibility 
of  Mr.  G.  F.  Mangum  of  the  Project  Management  Directorate.  Others  who 
cooperated  in  the  work  and  in  the  preparation  of  this  report  are  Messrs. 

E.  C.  Morganweck,  C.  S.  Combs,  C.  1.  Ashmore,  Dr.  W.  D.  Stephens  and 
Mrs.  E.  J.  Grice. 

(U)  This  report  has  been  assigned  the  Thiokol  internal  number  63-67{ControI 
No.  C-67-63A). 

(U)  This  project  is  being  accomplished  as  a  part  of  the  Air  Force  program, 
the  overall  objective  of  which  is  to  tailor  the  burning  rate  of  a  propellant 
predictably  and  concrollably  to  any  desired  level  in  the  range  from  1  to  10 
inches  per  second.  Experiments  will  be  performed  to  synthesize  more  ef¬ 
ficient  burn  rate  catalysts  by  maximizing  guideline  properties  determined 
under  Contract  AF04(6l  1)- 1 1212  toward  development  of  an  ideal  catalyst. 
Following  catalyst  synthesis,  each  compound  showing  potential  will  be  sub¬ 
jected  to  comprehensive  decomposition  studies  and  combustion  mechanism 
evaluation  with  propellant  ingredients. 

(U)  This  report  contains  no  classified  information  extracted  from  other 
classified  documents. 


STATEMENT  OF  APPROVAL 

(U)  Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


W.  H.  EBELKE,  Col.,  USAF 

Chief,  Propellants  Division 

Air  Force  Rocket  Propulsion  Laboratory 
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CONFIDENTIAL  ABSTRACT 


(C)  Thiokol's  program  ,o  tailor  the  burning  rate  of  a  propellant  predictably 
and  controllably  to  any  desired  level  in  the  range  from  1  to  10  inches  per 
second  is  divided  into  three  phases:  Phase  I  -  Synthesis  of  Burning  Rate 
Catalysts,  Phase  II  -  Decomposition  Studies  and  Evaluation  of  Catalysts  and 
Phase  III  -  Decomposition  of  Advanced  Oxidizers,  Fuels,  and  Binders. 
Experiments  will  be  performed  to  synthesize  more  efficient  burn  rate  catalysts 
by  maximizing  already  determined  guideline  properties  toward  development 
of  an  ideal  catalyst  under  Phases  I  and  II.  Following  catalyst  synthesis,  each 
compound  showing  potential  will  be  subjected  to  comprehensive  decomposition 
studies  and  combustion  mechanism  evaluation  with  propellant  ingredients. 
Thermogravimetric  analyses  of  catalysts  mixed  with  HC  polymer  show  that 
none  of  the  catalysts  tested  appear  to  catalyze  the  decomposition  of  HC  polymer. 
Propellants  formulated  with  the  copper  (I)  complex,  the  copper  (II)  complex  and 
the  cobalt  (II)  complex  of  0-ferrocenyl-£-oxo-propionaldehyde  demonstrated 
burning  rates  higher  than  that  of  similar  propellant  containing  n-butyl  ferrocene 
[PLASTISCAT-IV(R)]  at  the  2  percent  level.  An  analysis  to  relate  subsurface 
heating  to  changes  in  ammonium  perchlorate  crystal  structure  was  conducted. 

It  was  shown  that  if  the  concept  of  a  constant  surface  temperature  is  valid  and 
if  the  extrapolated  value  of  thermal  diffusivity  is  accurate  to  within  10  percent, 
the  surface  temperature  of  burning  ammonium  perchlorate  is  440  ±  30°C. 
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SECTION  1 


INTRODUCTION 


(U)  The  objective  of  this  program  is  to  tailor  the  burning  rate  of  a  propella  '■ 
predictably  and  controllably,  to  any  desired  level  in  the  range  from  1  to  10 
inches  per  second.  Ballistic  and  mechanical  properties  of  propellants  studied 
will  be  maintained  at  the  state-of-the-art  standards  of  current  Minuteman 
propellant. 

(C)  Thiokol's  approach  to  obtain  the  program  objective  will  be  through  the 
continued  investigation  of  iron  compound  effects  on  the  aluminum-ammonium 
perchlorate-polybutadiene  binder  system.  It  is  also  designed  to  systematical¬ 
ly  evaluate  new  propellant  ingredients  as  to  their  effect  on  burning  rate  and 
the  related  effect  on  combustion  mechanism.  New  materials  which  will  be 
evaluated  are  hydroxyl-ammonium  perchlorate,  hydrazine  diperchlorate  and 
nitronium  perchlorate  oxidizers;  aluminum  hydride,  beryllium,  beryllium 
hydride  fuels;  P-BEP,  NFPA  polymers,  and  TVOPA  plasticizer.  The  data 
and  conclusions  reached  under  Contract  AF04(6 1 1)- 1 12 12  will  be  used  as  a 
base  line  guide  for  the  work  to  be  accomplished  under  this  program. 

(U)  The  planned  program  consists  of  three  major  areas;  synthesis,  evaluation, 
and  advanced  ingredient  studies.  It  is  designed  to  progress  in  a  logical  manner 
so  that  the  most  promising  materials  receive  the  more  extensive  evaluation  and 
the  less  desirable  materials  are  discarded  at  an  early  date.  The  three  phases 
of  the  program  are: 

(U)  Phase  1  -  Synthesis  of  Burning  Rate  Catalysts 

(U)  Burning  rate  catalysts  will  be  synthesized  for  use  in  the  AP/A1/PB  pro¬ 
pellant  system.  Prior  knowledge  as  a  propulsion  contractor  and  the  data  ob¬ 
tained  from  Contract  AF04(6 1 1)- 1 12 1 2  will  serve  as  base  line  guides.  As  a 
new  material  is  synthesized,  it  will  be  evaluated  in  Phase  11. 

(U)  Phase  11  -  Decomposition  Studies  and  Evaluation  of  Catalysts 

(U)  An  evaluation  of  the  compatibility  of  the  candidate  catalysts  with  other 
propellant  ingredients  will  be  accomplished  upon  completion  of  synthesis, 
elemental  characterization,  and  physical  property  determination  of  each 
candidate  material.  Data  obtained  in  this  phase  will  allow  recommendations 
to  be  made  relative  to  the  development  of  new  burning  rate  catalysts  and  the 
development  of  high  burning  rate  propellants  utilizing  the  improved  catalyst. 
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(U)  Phase  III  -  Decomposition  of  Advanced  Oxidizers,  Fuels,  and  Binders 

(U)  Phase  III  effort  will  be  directed  toward  obtaining  a  basic  fundamental 
knowledge  of  the  decomposition  of  advanced  fuels,  oxidizers,  and  binders. 
Laboratory  test  data  will  be  utilized  to  postulate  a  burning  mechanism  of 
the  advanced  ingredient  and  a  comparison  made  with  that  of  conventional 
propellant  ingredients. 

(U)  This  report  covers  work  performed  for  the  period  1  June  through 
31  August  1967  under  Contract  F  04b  1  1  - 67- C- 0034.  Effort  to  date  has  been 
concerned  solely  with  Phases  I  and  II,  which  are  being  conducted  concurrently. 
Phase  III  will  be  initiated  after  the  completion  of  Phases  I  and  IL 


(This  pa 
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SECTION  II 


EXPERIMENTAL  ACCOMPLISHMENTS 

(U)  Based  cn  the  results  of  effort  expended  under  Contract  AF04(6l  1)- 1 1212, 
the  experiments  under  this  program  will  be  directed  toward  the  synthesis  and 
characterization  of  more  efficient  burn  rate  catalysts. 

(U)  i.  Phase  I  -  Synthesis  of  Burning  Rate  Catalysts 

(U)  7.'he  guide1  ines  for  synthesis  of  more  efficient  burn  rate  catalysts  have 
been  derived  from  effort  conducted  under  the  above  contract  and  are  listed 
below. 


High  iron  content 
Wide  liquid  range 
Readily  oxidizablc 

Compatible  with  other  propellant  ingredients 
Maximum  fuel  content  (heat  release  on  oxidation) 

(U)  The  synthesis  of  newer  burn  rate  catalysts  will  be  directed  toward  maxi¬ 
mizing  the  aforementioned  properties  into  the  ideal  catalyst. 


(U)  Effort  performed  in  this  area  of  research  during  this  reporting  period 
is  presented  in  subsequent  paragraphs. 


(C)  Esterification  of  8  -fer rocenoy lpropionic  acid  with  methanol  gave 
methyl  S  -ferrocenoyi  propionate  in  80  percent  yield. 


Reaction:  O 

Fc-C-CH2CH2COOH  +  CH^OH  - > 

O 

FcC-CH2CH2  COCH, 


(C)  Reaction  of  hydroxymethyl  ferrocene  with  4-hydroxy-2-butanone 
using  acetic  acid  catalyst  gave  4-fcrrocenylmethoxy-2-butanone  in 
77  percent  yield. 
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Reaction: 


Fc-CH2OH 


hochzch2< 


!CH. 


HOAc 


FcCH2OCH2CH2 


(C)  Substitution  of  2-methoxyethanol  for  4-hydroxy- 2-butanone  gave 
l-methoxy-2- (ferrocenylmethoxy)-ethane  in  70  percent  yield. 

Reaction: 


fc-ch2oh  +  hoch2ch2och3 - > 

fcch2och2ch2och3 

(C  )  Use  of  methylthioethanol  gave  1 -methylthio- 2-ferrocenylmethoxy- 
ethane,  in  50  percent  yield. 

Reaction: 

Fc-CH  OH  +  HOCH  CH  SCH.  - > 

FcCH2OCH2CH2SCH2 

(U)  Several  solid  derivatives  of  dimethylaminomethyl  ferrocene  have  been 
synthesized.  The  perchlorate,  nitrate,  picrate,  iodide,  bromide  and  fluoride 
salts  have  been  prepared  either  through  neutralization  or  metathesis  reactions. 

(C)  Treatment  of  N,  N- dimethylaminomethyl  ferrocene  with  bis- dimethyl- 
aminomethane  gave  1,  l'-bis  (dimethylaminomethyl)  ferrocene  in  28  percent 
yield. 


Reaction: 

CH  - > 

2 


FcCH2N(CH3)  + 
2 


(CH  )  N 
2 


Fc 


2 


(U)  The  structural  assignments  of  the  above  catalyst  were  confirmed 
by  infrared  analysis. 
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(C)  The  reaction  between  hydroxymethyl  ferrocene  and  mercaptoacetic 
acid  gave  «  - (fer rocenylmethylthio)-acetic  acid  in  75  percent  yield. 


Elemental  analysis: 


Carbon 

(%) 


Calculated  53.81 

F  ound  53.91 


Melting  point:  124  -  125°C 


Reaction: 


Hydrogen 

<%) 

4.  86 

5.  27 


Fc-CH  OH  +  HS-CH-COOH - >  Fc-CH  SCH  COOH 

2  c  2  £ 


(C)  Esterification  of  «  - (Ferrocenylmethythio)-acetic  acid  with  methanol 
gave  methyl  -  (Fer  roeenylmethylthio)  acetate,  according  to  the  following: 

Reaction : 


FcCH  SCH  COOH  +  CH-OH - >  Fc  CH-,  SCH_  COOCH, 

2  2  3  ^2  3 

(C)  The  reaction  of  acetyl  ferrocene  with  ethyl  mercaptan  gave  acetyl 
ferrocene  ethylthioketal  in  51  percent  yield. 

Melting  point:  40  -  4l°C 


Reaction: 

O 

fl 

Fc-C-CH  +  2CH3CH  SH 


(C)  An  unsuccessful  attempt  was  made  to  prepare  a  ke to- ester -sulfide 
ferrocene  compound  by  reaching  «  -  (ferrocenylmethylthio)- propionic  acid 
with  4-hydroxy- 2-butanone. 


CH, 

HC1  |  3 

- >  Fc-C-S-CH2CH3 

s-ch2ch3 
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(C)  A  ferrocene  compound  containing  keto-ester  and  ether  groups  was 
successfully  prepared  in  30  percent  yield  by  reacting  6  -ferrocenylpropionic 
acid  with  2-methoxyethanol,  which  gave  2  -  me  thoxye  thyl  -  S  -ferrocenoyl 
propionate. 

Melting  point:  47  -  49°C 
Reaction: 


O 

M 

Fc-C-CH  CH  COOH  + 
2  2 


O 


O 


HOCH2  CH2  OCH3— >Fc-C-CH2  CH2-C-OCH2  CH2  OCH3 


{C)  The  reaction  of  hydroxymethyl  ferrocene  with  2- chloroethanol  gave 
2-chloroethoxymethyl  ferrocene  in  50  percent  yield . 


Elemental  Analysis: 

Carbon 

Hydrogen 

(%) 

(%) 

Calculated 

56.1 

5.4 

F  ound 

57.9 

5.42 

Reaction:  Fc-CH2OH 

+  HOCH2CH2-Cl- 

— >  Fc-CH2OCH2CH2C1 

The  reaction  of  hydroxymethyl  ferrocene  with  ethyl  mercaptan 
thylthiomethyl  ferrocene  in  40  percent  yield. 

Elemental  Analysis: 

Carbon 

(%) 

Hydrogen 

(%) 

Calculated 

60.0 

6.16 

Found 

60.22 

6.33 

Reaction:  Fc-CH2OH  +  CH3CH2SH - *  Fc-CH2-S-CK2  CH3 

(C)  The  reaction  of  hydroxymethyl  ferrocene  with  propionic  acid  gave 
ferrocenylmethylpropionate  in  35  percent  yield. 
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Elemental  Analysis: 

Carbon 

Hydrogen 

<%) 

(%) 

Calculated 

61.  8 

5.  88 

F  ound 

62.  2 

5.  63 

Melting  point:  55  - 

57°C 

Reaction:  O 

O 

Fc-CH2OH+  HO-C-CH2CH3 - »  Fc-CH2-0-C-CH2CH3 

(C)  The  reaction  of  hydroxymethyl  ferrocene  with  hydrogen  sulfide 
gave  fe rrocenylmethyl  mercaptan  in  90  percent  yield. 

Melting  point:  53  -  55°C 

H+ 

Reaction:  Fc-CH2OH  +  H'2S  - - »  Fc-Cri2SH 

(C)  The  reaction  of  fer rocenylmethylmercaptan  with  chloro-2- 
propanone  in  sodium  hydroxide  solution  gave  l-(ferrocenylmethylthio)- 
2-propanone  in  65  percent  yield. 

O  O 

II  l> 

Reaction:  Fc-CH2S  Na  +  C1-CH2-C-CH3 - »  Fc-CH2-S-CH2-C-CH 

(C)  The  reaction  of  fe rr ocenylmethylmercaptan  with  allyl  bromide  in 
sodium  hydroxide  solution  gave  allylthiomethyl  ferrocene  in  80  percent 
yield. 


Reaction: 

Fc-CH2SNa  +  Br-CHZCH  =  CH2 _ >  Fc-CH2S-CH2-CH  =  CH2 

(C)  The  reaction  of  hydroxymethyl  ferrocene  with  isopropyl  alcohol 
gave  isopropoxymethyl  ferrocene  in  70  percent  yield. 

Melting  point:  33  -  34°C 


Reaction: 


CH_ 

I  3 


Fc-CH2OH  +  HOCH-CH3 - *■ 


CH 

I  3 

Fc-CH2-0-CH-CH3 


(C)  An  unsuccessful  attempt  was  made  to  prepare  a  keto-este r- sulfide 
ferrocene  compound  by  reacting  S -ferrocenoylpropionic  acid  with  2-methyl- 
thioethanol.  The  previously  attempted  esterification  of  a- (ferrocenyl- 
methylthio}- acecic  acid  with  4-hydroxy- 2-butanone  was  also  unsuccessful. 
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(U)  Physical  constants  ard  elemental  analysis  for  previously  reported 
catalysts  were  determined  during  this  report  period  and  are  shown  below. 


(C)  4-Ferrocenylmethoxy-2-butanone 

O 

II 

(Fc-CH  OCH  CH  C-CH3) 


Elemental  analysis: 


Carbon 


(%) 


Hydrogen 

(%) 


Calculated 

63.  0 

6.4 

Found 

63.  03 

6.42 

Melting  point:  15  -  17°C 

(C) 

1  -  me  thoxy-2-ferrocenylmethoxye  thane 

(Fc-ch2  och2ch 

2och3) 

Elemental  analysis: 

Carbon 

Hydrogen 

(%) 

(%) 

Calculated 

61.4 

6.  57 

F  ound 

61.  59 

6.  71 

Melting  point:  20  -  22°C 

(C) 

The  melting  point  of  methoxym< 

sthyl  ferrocene  (Fc-CH? 

determined  to  be:  6  -  8°C. 

CH, 

1  j 

(C) 

Acetyl  ferrocene  ethylthioketal 

[Fc-C  (-S- 

CH  -CH,)  ] 
2  3  2 

Elemental  analysis: 

Carbon 

Hydrogen 

(%) 

(%) 

Calculated 

57.  5 

6.6 

Fovnd 

57.  77 

6.  77 
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(C) 


o 

II 


2- me thoxy-  6 -ferrocenoyl  propionate  (Fe-C-CH  CH  C-O-CH,  CH?OCHJ 

o  2  c  &  3 


Elemental  analysis: 


Carbon 

Hydrogen 

(%) 

(%) 

Calculated 

59.4 

5.  82 

F  ound 

59.  5 

5.84 

The  melting  point  of  allyloxy  methyl  ferrocene  (Fe-CH^  OCH^CH  =  CH  ) 
was  determined  to  be  -  1 1  to  -  1 0°C.  ^ 


(U)  2.  Phase  11  -  Decomposition  Studies  and  Catalyst  Evaluation 

(U)  An  evaluation  of  the  compatibility  of  the  candidate  catalysts  with  other 
propellant  ingredients  will  be  accomplished  upon  completion  of  synthesis, 
elemental  characterization  and  physical  property  determination  of  each 
candidate  catalyst.  Following  the  compatibility  studies,  the  actual  effective¬ 
ness  of  the  candidate  catalysts  will  be  measured  by  determining  the  propellant 
processing  characteristics  as  well  as  the  burning  rate. 

(U)  a.  Catalyst  Compatibility  Studies 

(C)  Nineteen  newly  synthesized  catalysts  were  evaluated  by  differential 
thermal  analysis  (DTA)  during  this  report  period.  These  catalysts  were: 

2- methoxyethyl  ferrocenyl  ether 
hydroxymethyl  ferrocene 
2-methylthioethyl  ferr ocenylmethyl  ether 

A  ferrocene  polymer  of  undetermined  structure  (AFB  1-Poly) 

dimethyl  ferr  ocenylmethyl  ammonium  nitrate 

2,  2'  bis(dimethylaminomethyl)  ferrocene  dipicrate 

dime  thy  laminomethyl  ferrocene  picrate 

methyl- 4-  oxo-4-  ferr  ocenylbuty  rate 

acetyl  ferrocene  ethyithioketal 

cobalt  (II)  complex  of  6  -ferrocenyl-  8  -oxo-propionaldehyde 
ethylthiomethyl  ferrocene 

dimethyl  ferrocenylmethyl  ammonium  fluoride 
1-ferrocenylmethoxy-  3-butanone 
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manganese  (11)  complex  of  8  -ferrocenyl-  8  -oxo-propionaldehyde 

acetonylthiome thyl  ferrocene 

isopropyl  ferrocenylmethyl  ether 

2,  2-bis  (acetyl)  ethyl  ferrocene 

2-methoxyethyl-  8  -fer rocenoylpropionate 

2-chloroethyl  ferrocenylmethyl  ether 


These  data,  which  are  presented  on  Figures  1  through  6,  show  that  2-methyl- 
thioethyl  ferrocenylmethyl  ether,  dimethyl  ferrocenylmethyl  ammonium  nitrate, 
the  cobalt  (II)  complex  of  8  -ferrocenyl-  8  -oxo-propionaldehyde,  and  the 
manganese  (II)  complex  of  8 -ferrocenyi-  8  -oxo-propionaldehyde  appear  to 
affect  significantly  the  thermal  decomposition  of  ammonium  perchlorate. 


(C)  Several  previously  synthesized  catalysts  were  mixed  with  HC  polymer 
and  subjected  to  thermal  gravimetric  analysis  (TGA).  About  0.01  gram  of 
each  catalyst  was  mixed  with  0.  19  gram  HC  polymer,  producing  a  five  percent 
(weight)  catalyst  mixture.  HC  polymer  without  any  added  catalyst  was  also 
subjected  to  TGA.  Temperatures  at  the  start  and  at  the  end  of  weight  loss  are 
given  in  Tables  1  and  II.  The  lower  temperatures  at  start  of  weight  loss  for 
the  catalyst/HC  polymer  mixtures  are  attributed  to  volatilization  of  the  catalysts. 
The  catalyst/HC  polymer  mixtures  indicated  by  (a)  in  the  tables  apparently 
cured  to  a  rubber  during  the  pyrolysis.  The  first  sample  of  bis(  a -fe r rocenyl- 
ethyl)  ether/HC  polymer  was  examined  microscopically  after  being  thrown 
from  the  sample  tube.  It  was  a  hard  rubber  containing  numerous  bubbles  of 
gas.  The  audible  "pop"  accompanied  by  a  sudden  weight  loss,  noted  for  this 
and  some  of  the  other  samples,  is  therefore  attributed  to  sudden  release  of 
decomposition  gases  by  the  cured  HC  polymer.  Apparently  catalysts  containing 
an  ether  linkage  and  some  other  structures  react  to  cure  the  HC  polymer  at 
temperatures  around  360°C.  None  of  the  compounds  tested  appeared  to  catalyze 
the  decomposition  of  HC  polymer.  This  is  evidenced  by  the  fact  that  the  tem¬ 
perature  at  end  of  weight  loss  was  nearly  the  same  in  all  of  the  catalyst/HC 
polymer  mixtures. 


(C)  Recently  synthesized  catalysts  have  been  tested  for  compatibility  with 
HC  polymer  and  MAPO.  Small  scale  mixes  were  prepared  and  examined 
under  the  microscope  at  ambient  temperature  and  resulting  discoloration, 
dissolving  of  the  catalyst,  evolution  of  gas  or  heat,  etc.,  was  noted.  These 
samples  were  then  placed  in  a  1 6 0  _  1°F  oven  for  24  hours,  after  which  they 
were  re-examined  microscopically.  Control  tests  of  HC  polymer  and  MAPO 
containing  no  catalysts  were  also  run.  Results  are  given  in  Table  Ill.  All  of 
the  catalysts  tested  appear  to  be  compatible  with  HC  polymer  and  MAPO  at 
ambient  temperature.  However,  under  the  l60°F  conditions  only  2-methyl- 
thioethyl  ferrocenylmethyl  ether,  2-methoxyethyl  ferruceny  line  thyl  ether,  and 
dimethyl  ferrocenylmethyl  ammonium  nitrate  were  compatible  with  both  HC 
polymer  and  MAPO. 
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TIIERMOGRAVIMETRIC  ANALYSIS  OF  CATALYSTS  MIXED  WITH  HC  POLYMER  (U) 
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FiCH^SCH.CHj  --  ethyUkiuirurUivI  ferrocene 

FcCH.N(CH.|  •  HF  -  dimethyl  ierroceaylmethvl 
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Figure  4.  (C)  Differential  Thermal  Analyses  of  ethylthiomethyl  ferrocene,  dimethyl  ferrocenyl 
methyl  ammonium  fluoride,  and  l-ferrocenylmethoxy-3-butanone  . 
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(U)  b.  Burning  Rate  Studies 

(C)  A  1,000  gram  batch  of  HC  propellant  without  burning  rate  catalyst  or 
curing  agent  was  prepared,  and  100  gram  batches  of  propellant  containing 
burn  rate  catalyst  and  curing  agent  were  prepared  from  the  master  batch. 
Strands  were  cut  from  the  cured  propellant  and  burning  rates  were  determined. 
Catalysts  evaluated  were: 

copper  (II)  complex  of  6 -ferrocenyl-  £  - oxo-propionaldehyde 

60:40  mixture  of  1,  I-di  (methoxymethyl)  ferrocene  and 
1- hydroxymethyl-  1  '-methoxymethyl  ferrocene 

copper  (1)  complex  of  £  -ferrocenyl-  8  -  oxo-propionaldehyde 
dimethylaminomethyl  ferrocene  picrate 
n-butylth\omethyl  ferrocene 

cobalt  (II)  complex  of  8-fe rrocenyl-  £  -oxo-propionaldehyde 
manganese  (II)  complex  of  6-ferrocenyl-  8  -oxo-propionaldehyde 

(C)  Addition  of  copper  (II)  complex  followed  by  MAPO  addition  and  15  minutes 
of  vacuum  mixing  resulted  in  a  viscosity  of  60  to  65  kilopoise  at  145°F.  The 
addition  of  the  above  ether  mixture  gave  a  viscosity  reduction  to  20  kilopoise. 
No  viscosity  reduction  was  experienced  with  addition  of  either  the  copper  (I) 
complex,  dimethylaminomethyl  ferrocene  picrate,  the  cobalt  (II)  complex  or 
the  manganese  (II)  complex.  Addition  of  n- butylthiomethyl  ferrocene  to  the 
propellant  mix  effected  a  viscosity  reduction  to  the  20  to  25  kilopoise  range. 

(C)  Propellant  made  with  the  copper  (II)  complex  gave  a  very  soft  cure  after 
approximately  one  month.  A  good  cure  was  exhibited  by  the  propellant  made 
with  the  ether  mixture  after  four- and- one /half- days.  Propellant  made  with 
the  copper  (I)  complex  gave  a  soft  cure  after  approximately  one  month.  The 
degree  of  cure  exhibited  was  similar  to  that  for  the  copper  (II)  complex  of  the 
same  compound.  Dimethylaminomethyl  ferrocene  picrate  gave  propellant  that 
was  hard  and  crumbly.  A  good  propellant  cure  was  produced  by  n- butylthio¬ 
methyl  ferrocene.  Propellant  made  with  the  cobalt  (II)  complex  cured  in 
about  three  days  and  was  hard  and  crumbly.  It  is  probable  that  the  incom¬ 
patibility  of  the  cooalt  (11)  complex  with  HC  polymer  at  elevated  temperatures 
was  the  cause  of  this  "hard"  cure.  Propellant  made  with  the  manganese  (II) 
complex  gave  a  good  propellant  cure. 
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(C)  A  comparison  of  the  burning  rate  of  propellant  formulated  with  these 
complexes  and  with  n-butyl  ferrocene  [ PLASTISCAT-IV^ )  ]  (Figures  7 
through  11  )  show  the  following: 


Copper  (II)  complex 
60:40  ether  mixture 

Copper  (I)  complex 

Dimethylaminomethyl  ferrocene 
picrate 

n-butylthiomethyl  ferrocene 
Cobalt  (II)  complex 

Manganese  (II)  complex 


Higher  over  entire  pressure  range 

Slightly  higher  over  entire  pressure 
range 

Considerably  higher  throughout  pressure 
range  of  400  to  1200  psig 

Considerably  below 


Considerably  below 

Somewhat  higher  throughout  pressure 
range  of  40f  tc  1200  psig 

Somewhat  lower  throughout  pressure 
range  of  400  to  1200  psig 


These  data  also  show  that  burning  rates  for  the  copper  (I)  complex  are  higher 
than  the  burning  rates  for  the  copper  (II)  complex  of  the  same  parent  compound. 


(U)  c.  Analysis  to  Relate  Subsurface  Heating  to  Changes  in 

Ammonium  Perchlorate  Crystal  Structure 

(C)  An  analysis  to  relate  subsurface  heating  to  changes  in  ammonium  per- 
chorate  crystal  structure  was  conducted  and  is  included  as  Appendix  A  to  this 
report.  It  was  concluded  from  this  analysis  that,  if  the  concept  of  a.  constant 
surface  temperature  is  valid  and  if  the  extrapolated  value  of  thermal  diffusivity 
is  accurate  to  within  10  percent,  the  surface  temperature  of  burning  ammonium 
perchlorate  is  440  1  30°C.  In  addition,  calculations  made  using  very  high 
heat  generation  rates  indicate  that  subsurface  reactions  are  not  an  important 
factor  in  controlling  the  burning  rate  of  ammonium  perchlorate  at  pressures 
above  500  psi. 
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SECTION  m 

CONCLUSIONS 


(C)  The  following  conclusions  may  be  drawn  from  the  work  conducted 
during  this  reporting  period: 

Thermogravimetric  analyses  (TGA)  of  catalysts  mixed  with  K  C 
polymer  showed  that  none  of  the  compounds  tested  was  effective 
in  catalyzing  the  decomposition  of  HC  polymer. 

All  of  the  catalysts  tested  appear  to  be  compatible  with  HC 
polymer  and  MAPO  at  ambient  temperature;  however,  under 
the  160°F  conditions,  only  2-methylthioethyl  ferrocenylmethyl 
ether,  2-methoxyethyl  ferrocenylmethyl  ether,  and  dimethyl 
ferrocenylmethyl  ammonium  nitrate  were  compatible  with 
both  HC  polymer  and  MAPO  . 

Propellants  formulated  with  the  copper  (I)  complex,  the 
copper  (II)  complex  and  the  cobalt  (n)  complex  of 
0-ferrocenvl-  0-oxo- propionaldehyde  demonstrated  burning 
rates  higher  than  that  of  similar  propellant  containing  n-butyl 
ferrocene  [PLASTISCAT-IV(R)]  at  the  2  percent  level. 

The  burning  rate  of  propellant  formulated  with  the  copper  (I) 
complex  was  higher  than  the  burning  rate  for  the  copper  (II) 
complex  of  the  same  parent  material. 

The  sulfur  containing  ferrocene  derivatives  synthesized 
were  for  the  most  part  solids  and  were  not  effective 
as  burn  rate  modifiers. 

The  salts  of  dimethylaminomethyl  ferrocene  do  not  contain 
enough  iron  to  be  effective  solid  bum  rate  catalysts  at  levels 
of  1  to  5  percent. 

An  analysis  to  relate  subsurface  heating  to  changes  in 
ammonium  perchlorate  crystal  structure  was  conducted.  It 
was  concluded  from  this  analysis  that,  if  the  concept  of  a 
constant  surface  temperature  is  valid  and  if  the  extrapolated 
value  of  thermal  diffusivity  is  accurate  to  within  10  percent, 
the  surface  temperature  of  burning  ammonium  perchlorate 
is  440  1  3  0° C.  In  addition,  calculations  made  using  high  heat 
generation  rates  indicate  that  subsurface  reactions  are  not  an 
important  factor  in  controlling  the  burning  rate  of  ammonium 
perchlorate  at  pressures  above  5  00  psi. 
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APPENDIX  A 


(Unclassified) 

ANALYSIS  TO  RELATE  SUBSURFACE  HEATING  TO  CHANGES 
IN  AMMONIUM  PERCHLORATE  CRYSTAL  STRUCTURE 


(U)  1.  Introduction 

(U)  An  analysis  was  conducted  to  relate  surface  temperature,  subsurface 
reactions,  and  crystal  structure  of  burning  ammonium  perchlorate  (AP)  using 
a  more  complete  mathematical  model  to  analyze  further  the  data  presented  in 
a  series  of  papers  by  Beckstead,  Hightower,  and  Price^"^.  The  results  of 
this  analysis  have  direct  application  to  the  development  and  application  of 
burning  rate  models  to  explain  observed  propellant  phenomena.  Since  the 
burning  rate  of  heterogeneous  propellants  containing  AP  is  controlled  by  the 
AP  particle  size  and  percentage,  it  is. not  sufficient  to  treat  propellant 
burning  in  the  mean  but  rather  the  burning  of  individual  ingredients  must  be 
considered.  To  account  accurately  for  the  heat  transferred  to  the  burning 
AP  surface,  the  contribution  of  subsurface  heating  and  the  value  of  surface 
temperature  must  be  known.  Since  AP  decomposes  at  above  250°C  and  the 
surface  temperature  of  burning  composite  propellants  containing  AP  is 
above  400°C,  considerable  difficulty  is  being  encountered  in  measuring  its 
thermophysical  characteristics  in  the  temperature  range  where  the  burning 
rate  controlling  reactions  occur,  i.  e.  ,  above  250°C.  Several  carefully 

ii  ^ 

controlled  optical  measurements  ’  and  embedded  micr o-thermocouple 
measurements  ^ ^  Qf  surface  temperature  have  been  conducted.  These 
experimental  techniques  were  recently  reviewed  by  Friedman^.  In  general, 
the  optical  teenniques  are  not  effective  at  rocket  motor  operating  pressures 
and  micr  o-thermocouple  s  are  too  large  to  follow  the  temperature  rises  in 
the  propeilant  ingredients.  In  addition,  since  the  propellant  is  a  heterogeneous 
substance,  difficulty  is  encountered  in  discriminating  between  mean  and  local 
effects. 

(U,  The  basis  of  this  analysis  is  the  establishment  of  the  heat  balance  to 
correlate  experimental  observations  of  the  crystallographic  phase  transition 
from  an  orthorhombic  to  a  cubic  structure  in  burning  AP  crystals. 


(U)  2.  Background 

IU)  Several  investigators  have  recognized  that  since  the  crystalline  trans¬ 
formation  from  orthorhombic  to  cubic  occurs  approximately  midway  between 
the  ambient  propellant  temperature  and  the  burning  surface  temperature, 
it  cc.’Jd  be  used  to  deduce  surface  temperatures  of  burning  propellants. 
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McGurk  examined  optically  very  thin  microtomed  sections  of  AP  composite 
propellants  before  and  after  extinguishment.  He  noted  an  optical  discontinuity 
below  the  extinguished  surface  of  individual  AP  particles  which  he  postulated 
to  be  the  interface  between  the  orthorhombic  and  cubic  phases  prior  to 
propellant  extinguishment.  Of  course,  after  extinguishment,  the  cubic 
structure  reverts  back  to  the  orthorhombic  structure.  In  addition,  McGurk 
suggested  that  the  100°C  isotherm  could  be  bracketed  by  examining  the  effect 
of  surface  heating  on  voids  containing  a  saturated  water  and  ammonium 
perrnlorate  solution.  However,  McGurk  did  not  make  quantitative  deter¬ 
minations  of  surface  temperature. 

(U)  Selzer  photographed  thin  slices  of  burning  composite  propellant  and 
established  the  cubic  phase  thickness  of  AP  from  patterns  produced  by 
polarized  light.  Although  Selzer 's  technique  produced  a  wealth  of  data,  his 
attempt  to  deduce  a  value  of  surface  temperature  was  premature  because 
of  experimental  and  analytical  difficulties.  Since  his  optics  were  penetrating 
a  relatively  thick  section  of  propellant,  he  was  not  measuring  the  cubic  phase 
thickness  at  a  discrete  point.  In  addition, as  pointed  out  by  Beckstead  and 
Hightower^,  Selzer's  analytical  treatment  to  deduce  surface  temperature 
was  not  sufficiently  comprehensive. 

(U)  Because  of  the  transient  nature  of  the  burning  of  individual  AP  particles, 
it  is  doubtful  the  data  of  McGurk  and  Selzer  can  be  used  to  make  accurate 
determinations  of  the  heat  balance  that  occurs  in  the  solid  phase.  As  part  of 
the  work^  being  conducted  to  develop  a  more  comprehensive  burning  rate 
model,  equations  were  solved  that  approximate  the  transient  two-dimensional 
(the  third  dimension  is  accounted  for  by  spherical  symmetry)  heat  conduction 
that  occurs  as  individual  AP  particles  are  consumed.  Figure  1  shews  ana¬ 
lytical  results  for  three  spherical  AP  particles  being  burned  at  700  psi. 
Similar  transients  have  been  observed  in  micro-window  bomb  photographs. 
The  calcuations  reveal  that  because  the  AP  particles  are  preheated  as  the 
adjacent  binder  layer  is  consumed,  a  burning  rate  overshoot  occurs.  This 
initial  overshoot  decays  as  the  preheated  portion  of  the.  AP  is  consumed. 

As  the  spherical  AP  particle  burns,  the  ratio  of  exposed  surface  to  immedi¬ 
ate  subsurface  volume  continuously  varies.  This  coupled  with  the  normal 
induction  time  to  establish  thermal  profiles  results  in  continuously  varying 
burning  rates.  Just  prior  to  the  burnout  of  the  AP  particle,  heat  is  trans¬ 
ferred  to  the  binder  layer  that  separates  adjacent  AP  particles.  Since 
binder  has  a  thermal  conductivity  approximately  50%  lower  than  the  thermal 
conductivity  of  AP,  moi  e  heat  is  retained  in  the  AP  crystal  causing  the 
rapid  increase  in  burning  rate. 

(U)  The  excellent  experiments  recently  reported  by  Beckstead  and  Hightower 

'  w>  rvinnu  nf  fKo  cVmffmrviinrrc  nf  fUo  a  i  r  i  nil  o  o  c?  V*  i  y  o  A 
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very  pure,  single  AP  crystals  in  a  nitrogen  atmosphere.  The  crystals  were 
quenched  by  rapid  depressurization.  The  cubic  phase  thickness  was  estab¬ 
lished  by  examining  many  cleaved  crystals  under  polarized  light  (Figure  2). 
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FIGURE  A-l.  BURNING  TRANSIENTS  OF  SPHERICAL  AMMONIUM 
PERCHLORATE  WITH  PARTICLE  DIAMETER  AS  A 
PARAMETER 
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CUBIC  PHASE  THICKNESS,  x  microns 
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RECIPROCAL  OF  CUBIC  PHASE  THICKNESS, 

1/x  ,  microns 
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FIGURE  A- 2.  THE  PHASE  TRANSITION  THICKNESS  AS  A 

FUNCTION  OF  BURNING  RATE  AS  MEASURED 
BY  BECKSTEAD  AND  HIGHTOWER2 
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A  very  important  conclusion  was  that  the  product  of  AP  burning  rate,  r, 
times  cubic  phase  thickness,  x  ^  p ,  was  nearly  constant.  In  addition,  the 
circles  on  Figure  2  are  indicative  of  the  variations  in  measuring  a  large 
number  of  samples^,  hence  the  scatter  observed  by  Selzer  did  not  occur 
in  the  single  crystal  experiments.  Beckstead  and  Hightower ^  solved  a 
boundary  value  problem  for  the  heat  balance  at  the  c  ibic-ortho- 
rhombic  interface  and  arrived  at  the  important  conclusion  that  the  tempera¬ 
ture  cf  the  solid  phase  surface  is  nearly  constant.  They  postulated  that  the 
higher  surface  temperatures  measured  optically’’  were  measurements  of  the 
surface  of  a  liquid  layer  which  Hightower  had  previously  reported  to  exist^. 

(U)  3.  Analysis 

(U)  As  shown  on  Figure  3,  the  decomposition  of  a  burning  AP  crystal  can  be 

considered  to  occur  in  several  phases:  1)  relatively  gradual  heating  of  the 

orthorhombic  phase  (<  1  0^  °C /sec  for  r  =  0.  8  cm /sec),  2)  transition  from 

orthorhombic  structure  to  cubic  structure,  3)  relatively  rapid  heating  of  the 

cubic  phase  (>10^  °C/sec),  4)  the  possible  formation  of  a  liquid  layer, 

5)  sublimation  or  gasification  to  AP  vapors,  and  6)  gas  phase  decomposition 

of  AP  vapors.  Since  for  steady-state  conditions  the  surface  of  the  solid 

phase  is  regressing  at  a  rate  equal  to  the  propagation  rate  of  the  transition 

zone,  the  thickness  of  the  cubic  phase,  x  ,  is  constant.  The  transition 

t 

from  cubic  to  orthorhombic  is  an  endothermic  process  which  requires  19.  5 
cal/gram  at  slow  heating  rates  (less  than  100°C  per  minute).  Since  this 
heating  rate  is  very  much  Jess  than  the  10^  °C/sec  that  occurs  during  normal 
propellant  burning,  it  is  not  entirely  certain  that  the  same  endothermic 
reaction  occurs  in  a  propellant  burning  at  pressures  above  500  psi.  Indeed, 
one  of  the  stated  objectives  of  this  analysis  is  to  evaluate  the  extent  to  which 
AP  subsurface  reactions  are  a  factor  in  controlling  propellant  burning  rate. 
By  selecting  the  AP  solid  phase  boundary  as  a  reference  point,  the  effect  of 
the  postulated  liquid  phase  is  not  a  factor  in  this  analysis. 

(U)  While  the  relative  effects  can  be  studied  using  approximate  values  for 
the  thermal  properties,  the  objective  of  this  analysis  is  to  obtain  quantitative 
results  by  taking  full  advantage  of  the  best  available  thermal  property 
determinations.  Since  orthorhombic  phase  thermal  properties  are  relatively 
well  established  compared  to  those  ot  the  cubic  phase,  the  thermal  property 
discussion  will  be  limited  to  the  cubic  phase.  The  cubic  phase  density  of 
1.  76  gm/cc  calculated  by  Marlcowitz  ^  was  selected.  However,  hot  stage 
photographs  of  single  AP  crystals  above  260°C  taken  by  Flanigan  have  shown 
that  the  crystal  growth  does  not  support  the  change  in  density  predicted  by 
Markowitz  .  The  specific  heat  values  were  taken  from  the  JANAF  data^. 
The  discontinuity  in  specific  heat  that  occurs  at  243°C  is  consistent  with  the 
predicted  change  in  density.  The  orthorhombic  phase  thermal  conductivity 
of  Rosser,  lnami,  and  Wisel6  provides  a  basis  for  determing  thermal 
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FIGURE  A-3.  PHYSICAL  SITUATION  CONSIDERED  BY 
MATHEMATICAL  SOLUTION 
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conductivity  of  the  cubic  phase  by  extrapolation.  The  problem  is  to  decide 
how  to  extrapolate  the  data  to  the  cubic  phase.  The  data  points  on  Figures 
4  and  5  are  Rosser's  data  corrected  for  porosity  in  the  manner  suggested  in 
Reference  16.  Beckstead  and  Hightower^  consider  the  data  point  at  215°C 
to  be  suspect.  If  it  is  ignored,  the  data  can  be  correlated  as  a  function  of 
reciprocal  absolute  temperature.  The  thermal  conductivities  of  other  crystals 
have  been  correlated  over  a  wide  range  of  temperatures  using  reciprocal 
temperature.  Thus,  by  extrapolating  the  data  a  value  of  thermal  conductivity 
at  400  C  of  0.00084  cal/ sec-cm- °C  is  obtained  from  the  linear  correlation 
{Figure  4)  and  a  value  of  0.  00097  cal/sec-cm-°C  is  obtained  from  the  recip¬ 
rocal  temperature  correlation  {Figure  5).  A  -10%  error  in  cubic  phase 
thermal  conductivity  will  increase  calculated  surface  temperature  by  approxi¬ 
mately  30 °C.  The  arguments  put  forth  by  Beckstead  and  Hightower  for 
extrapolation  based  on  reciprocal  temperature  appear  reasonable  and  will  be 
adopted  for  this  analysis. 

{U)  a.  Conduction  in  Ammonium  Perchlorate 

{U)  Using  the  physical  situation  described  above,  the  equations  governing 
conduction  within  a  moving  homogeneous,  isotropic,  reacting  medium  is 

P  C  dT/dt  -  a  [k  {  d  T/a  x)  j  /  a  X  +  r  3  T/a  x  +  q'  {1) 

The  sign  convention  for  burning  AP  is:  burning  rate,  r,  and  latent  heat  of 
crystallographic  phase  transition  (for  endothermic  reaction),  X  ,  less  than 
zero.  During  steady- state  burning  of  the  AP  crystal  the  boundary  conditions 
are: 


At  x  =  x  T  =  T 

t  r  t  r 

-k  dT/dx-L-k2  dT/dx]=rp2X  (2) 

After  extinguishment  the  boundary  conditions  become: 


At  x  =  0: 

T  < 

T 

s 

At  x  -  x 

t  r 

T  = 

T  - 
t  r 

r  4-  0 

T  < 

T  > 
t  r 

r  =  0 

The  above  equations  were  solved  using  explicit  finite  difference  techniques. 

t>  -  _  _  . .  n  _  _  r  *.i-  _ _ _ . _ -  *  n  a.  i  ~  j:.*. - -  : _ _  ~  ...  a.  .  4.1.  *  j. _ ~  ~  .  j 

l_>  c  e  a  u  &  c  ui  tiic  vex  y  oxuau  time  anu  Giotauuc  mticiuciito  ciiat  were  uocu, 

considerable  attention  had  to  be  given  to  developing  an  efficient  computational 
procedure.  The  steady-state  profiles  were  established  by  relaxation; 
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FIGURE  A-4.  THERMAL  CONDUCTIVITY  AND  SPECIFIC 
HEAT  OF  AMMONIUM  PERCHLORATE 
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THERMAL  CONDUCTIVITY,  k  x  10  ,  cal/sec  cm 
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TEMPERATURE,  T,  °C 


RECIPROCAL  TEMPERATURE,  1/T“K 

FIGURE  A-5.  RELATIONSHIP  BETWEEN  THERMAL 
CONDUCTIVITY  AND  TEMPERATURE 
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the  temperature  time  history  after  extinguishment  was  calculated  using  a 
forward  integration. 

(U)  b.  Steady-State  Closed  Form  Solutions 

(U)  Assuming  steady-state  conditions,  constant  thermal  properties,  and  no 
heat  generation!  Equation  1  reduces  to 

0  =  (k/p  c)  (d^T  /d  x^)  +  r  T /d  x  (3) 


Solving  Equation  3  with  the  steady-state  boundary  conditions,  yields 


’  Ttr  +  [c*p  (-rxti/ai) 


-  1]  t  <P  ,  c 


2  2 


/P^j) 


t  r 


(P2  C2/P1C1>  TC  '  f)2X/plcl  ]  <4) 

This  equation  differs  from  the  equation  reported  by  Beckstead  and  Hightower 
because  of  the  manner  in  which  they  considered  the  boundary  condition  at  the 
cubic -orthorhombic  interface.  They  wrote  boundary  condition  corresponding 
to  Equation  3  as 


-k  dT/'dx  -  [- 


k2  dT/dx] 


+  r  T 


t  r 


(p2C2 


Vl>  = 


r  P, 


(5) 


The  last  term  on  the  left  hand  side  of  the  equation  was  included  to  account 
for  the  change  in  heat  content  across  the  cubic-orthor hcmbic  interface.  This 
term  is  considered  to  be  redundant  because  the  effect  is  accounted  for  in  the 
value  of  latent  heat,  X  .  The  equation  for  steady- state  condition  derived  by 
Beckstead  and  Hightower  is 

Ts  =  Ttr  +L.xp(-ixtr/ai)  -UtTtr-  (P2p2/P,p1)T0-P2^/P1c1  ]  (61 
(U)  c.  Effect  of  Physical  Property  Accuracy 

(in  To  evaluate  the  effect  of  physical  property  accuracy  on  calculated  surface 
temperature,  the  variables  in  Equation  4  were  perturbed.  The  results  of 
Figure  6  were  obtained  using  the  indicated  properties  to  define  the  datum  case. 
The  results  indicate  that  the  property  which  is  least  understood,  the  cubic 
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SURFACE  TEMPERATURE,  T  ,  °C 

s 

FIGURE  A-6 .  EFFECT  OF  PHYSICAL  PROPERTY  ACCURACY 
ON  SURFACE  TEMPERATURE 
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phase  thermal  diffusivity,  has  the  greatest  effect  on  the  accuracy  of  the 
calculated  surface  temperature  value.  In  addition,  the  values  of  cj  and 
Pj  are  also  uncertain.  Since  Cj  andp^  are  always  used  as  a  product,  p^c^, 
it  is  felt  that  the  value  of  the  product  is  known  to  within  5%.  The 
accuracy  of  the  other  physical  property  values  used  in  the  equation  are 
well  known  relative  to  the  accuracy  of  a  . 

(U)  3.  Results 

(U)  A  series  of  numerical  computations  were  performed  in  which  the  values 
of  the  solid  phase  surface  temperature  and  burning  i  ate  were  varied.  The 
first  calculations  considered  that  the  only  subsurface  reaction  was  during  the 
transition  from  orthorhombic  to  cubic  crystal  structure;  the  thermal  physical 
properties  of  Figure  5  were  used.  Figure  7  indicates  the  relative  importance 
of  the  transients  after  strand  extinguishment.  For  a  surface  regression  of 
0.  6  cm/sec,  for  the  lower  value  of  surface  temperature  (410°C)  approximately 
0.  5  millisecond  is  required  for  the  regression  of  the  cubic  layer  to  stop.  For 
the  higher  surface  temperature,  460°C,  approximately  0.7  millisecond  is 
required.  As  the  regression  rate  increases  from  0.  6  cm/sec  to  1 .  Z  cm/sec, 
the  increase  in  the  cubic  phase  thickness  that  occurs  after  extinguishment 
varies  from  1 1  to  10%,  respectively. 

(U)  Figure  8  shows  the  effect  of  surface  temperature  on  the  thickness  of  the 
cubic  phase.  When  the  surface  temperature  is  related  to  the  measured  cubic 
phase  thickness,  a  value  of  470°C  is  obtained  if  the  transients  after  extinguish¬ 
ment  are  ignored,  However,  if  the  measured  cubic  phase  thickness  is  related 
to  the  predicted  cubic  phase  thickness  after  cooldown,  a  decreased  predicted 
surface  temperature  of  44Z°C  is  obtained. 

(U)  The  difference  between  the  exact  and  numerical  solution  foi  the  steady- 
state  condition  is  less  than  1%.  Using  the  variable  thermal  properties  of 
Figures  4  and  5  results  in  slightly  higher  values  for  the  steady- state  cubic 
phase  thickness  and  the  same  value  for  the  cubic  phase  thickness  after 
cooldown  'Figure  8).  This  agreement  is  expected  because  for  steady- 
state  the  net  heat  balance  can  be  accounted  for  exactly  by  the  mean  thermal 
properties;  and  for  the  transient  condition  the  variation  in  temperature  is 
small  enough  that  the  properties  evaluated  at  mean  temperatures  accurately 
represent  the  heat  balance. 

(U)  Figure  9  shows  cubic  phase  thickness  over  a  range  of  burning  rates. 

The  circle^  are  deduced  values  of  surface  temperature  based  on  the  measured 
cubic  phase  thicknesses. 
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FIGURE  A 


7.  TRANSIENTS  AFTER  EXTINGUISHMENT 
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STEADY  STATE  ,  r  -  .  8  cm/sec 
AFTER  COOL  DOWN 


SURFACE  TEMPERATURE,  T  ,  °C 
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FIGURE  A-8.  CUBIC  PHASE  THICKNESS  AS  A  FUNCTION 
OF  SURFACE  TEMPERATURE 
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CUBIC  PHASE  THICKNESS  AND  SURFACE 
TEMPERATURE  OVER  A  RANGE  OF 
SURFACE  TEMPERATURE 
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(U)  a.  Subsurface  Reactions 

(U)  Figure  10  is  a  cross-plot  ol'  the  surface  temperatures  of  Figure  8  and 
shows  that  surface  temperature  increases  slightly  with  increasing  burning 
rate.  The  slight  increase  in  T  provides  a  basis  for  discussing  the  rela¬ 
tive  eflects  of  subsurface  reactions.  It  will  be  assumed  that  the  surface 
temperature  of  the  solid  phase  is  constant  and  that  the  increase  in  calcu¬ 
lated  surface  temperature  with  increasing  rate  can  be  explained  by  the 
experimental  error  in  the  measurement  of  .  The  assumption  of  a 
constant  surface  temperature  is  not  inconsistent  with  data  in  light  of  the 
observations  made  by  Beckstead  and  Hightower^  that  surface  temperatures 
measured  optically  is  the  surface  temperatures  at  the  surface  of  a  liquid  phase. 

(U)  Two  conditions  may  be  considered:  1)  a  net  subsurface  endothermic 
reaction,  and  2)  a  net  subsurface  exothermic  reaction.  By  net  reaction 
the  latent  heat  of  phase  change  between  the  orthorhombic  and  cubic  crystal 
structure  is  included  with  any  subsurface  reaction  that  may  occur  in  the  cubic 
phase.  Figure  11a  shows  that  diminishing  the  net  endothermic  reaction  (or 
regression  toward  a  greater  exothermic  reaction)  the  actual  surface  tempera¬ 
ture  will  decrease.  Also  if  there  were  a  subsurface  reaction,  its  fraction  of 
depletion  would  tend  to  diminish  with  increasing  rate  as  indicated  in  Figure  lib. 

(U)  Now  we  will  consider  the  situation;  1)  during  the  actual  burning  of  AP 

there  is  a  net  solid  phase  endothermic  reaction,  and  it  is  ignored  in  the 

calculations,  and  2)  the  measured  values  of  x  were  used  to  deduce  surface 

t  r 

temperature.  The  trends  shown  in  Figures  11a  and  lib  indicate  that  as  a 
result  of  the  above  situation,  the  error  in  the  deduced  value  of  surface  temp¬ 
erature  wculd  be  in  the  direction  shown  on  Figure  10.  If  a  similar  argument 
is  attempted  for  a  net  exothermic  reaction,  the  reverse  trend  of  that  indi¬ 
cated  by  Figure  10  is  required.  Thus,  the  surface  temperature  variations 
of  Figure  10  are  not  great  enough  to  conclusively  support  a  net  endotherm, 
they  certainly  do  not  support  the  existence  of  an  exothermic  reaction.  In 
any  event,  the  subsurface  reaction  is  not  a  primary  effect,  for,  if  it  were, 
greater  discrepancy  in  calculated  (or  deduced)  values  of  surface  temperature 
would  occur. 

(U)  The  calculations  of  the  cubic  phase  thickness  were  repeated  with  the 
subsurface  heat  generation  expressed  by  the  first  order  reaction  equation 

qi  =  Z  Q  f  (t)  exp  (-E/RT) 

The  heat  evolution  in  the  solid  phase  was  assumed  to  be  200  cal/gm  and  the 
value  of  the  frequency  factor  was  varied  from  10^  to  10^®,  Z  equal  to  10^ 
is  considerably  above  the  highest  values^7  that  have  been  assigned  to  it. 

Even  for  the  high  values  of  frequency  factor,  less  than  1%  change  in  surface 
temperature  was  calculated.  While  using  these  external  assumptions,  a 
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subsurface  heat  release  within  the  first  two  microns  was  calculated; 
the  heat  release  was  not  of  sufficient  magnitude  to  affect  burning  rate. 

Again,  these  calculations  confirm  that  even  though  the  solid  phase 
temperatures  are  of  the  proper  magnitude  for  subsurface  reactions 
to  occur,  the  stay  time  is  too  short  for  the  reaction  to  the  extent  that 
it  affects  burning  rate. 

(U)  On  the  basis  of  kinetic  data  obtained  at  very  low  heating  rates  {less  than 
100°C/min)  Feinauer,  Waesche,  and  Wenograd^  have  postulated  that  the 
condensed  phase  reactions  play  an  important  role  in  the  control  of  composite 
solid  propellant  burning  rates.  The  results  of  this  analysis  conclusively 
demonstrate  that  since  the  time  available  for  subsurface  reactions  to  occur 
in  practical  propellants  is  several  orders  of  magnitude  less  than  the  heating 
periods  used  in  the  laboratory  apparatus  employed  to  obtain  the  kinetic  data, 
the  condensed  phase  reactions  are  not  of  sufficient  magnitude  to  control 
propellant  burning  rate. 

(U)  4.  Conclusions 

(C)  If  the  concept  of  a  constant  surtace  temperature  is  valid  and  if  the  extra¬ 
polated  value  of  thermal  diffusivity  is  accurate  to  within  10%,  the  surface 
temperature  of  burning  AP  is  440  ±  30 °C.  In  addition,  calculations  made 
using  very  high  heat  generation  rates  indicate  that  subsurface  reactions  are 
not  an  important  factor  in  controlling  the  burning  rate  of  AP  at  pressures 
above  500  psi. 
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Latin  Symbols 


Definition 


E 

f 

k 

'q 

Q 


R 


specific  heat 

activation  energy  for  pyrolysis  law 
fraction  of  reactant  remaining  unconsumed 
thermal  conductivity 

energy  per  unit  volume  from  chemical 
reactions 

heat  evolution  accompanying  a  chemical 
or  physical  change 

burning  rate 

gas  constant 


T 

t 


Greek  Symbols 

a 

\ 

P 


temperature 

time 

normal  distance  from  the  surface 
pre-exponential  factor 

thermal  diffusivity 

latent  heat  of  crystal  phase  transition 
density 


Subscripts 

s 

tr 


solid  phase  surface 
cubic-orthorhombic  interface 


Unit 

cal/gm°F 

cal/mole 

cal/cm  sec°C 

O 

cal/cm  sec 
cal/gm 

cm/  sec 

1. 9867 

cal/mole°K 

°K 

sec 

cm 

2  , 

cm  /sec 
cal  /  gm 
gm/cru^ 
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